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Abstract 

Using the overlap-Dirac operator proposed by Neuberger, we have computed in 
lattice QCD the one-loop renormalization factors of ten operators which measure 
the lowest two moments of unpolarized and polarized non-singlet quark distribu- 
tions. These factors are necessary to extract physical numbers from Monte Carlo 
simulations made with overlap fermions. 

An exact chiral symmetry is maintained in all our results, and the renormalization 
constants of corresponding unpolarized and polarized operators which differ by a 75 
matrix have the same value. We have considered two lattice representations for each 
continuum operator. The computations have been carried out using the symbolic 
language FORM, in a general covariant gauge. In some simple cases they have also 
been checked by hand. 
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1 Introduction 



In recent years remarkable progress has been made towards formulations of 
lattice fermions which have no doublers and possess an exact chiral symmetry 
without giving up desirable features like flavor symmetry, locality, unitarity 
or gauge- invariance. The Ginsparg- Wilson relation [1] 

l5 D + D l5 = a-D l5 D (1) 
P 



has been recognized [2] as the fundamental condition in this context, and 
a Dirac operator D which satisfies this relation can indeed describe chiral 
fermions on the lattice. One of the possible solutions of the Ginsparg- Wilson 
relation has been found by Neuberger [3] starting from the overlap formal- 
ism [4]. In the massless case the overlap-Dirac operator is 



D 



N 



-P 



1 + 



X - 
y/XXX- 



(2) 



where 



X = D w - -p (3) 

a 



in terms of the usual Wilson-Dirac operator 

1 



D 



w 



V ij,ip(x) = - U(x, n)i){x + afi) - tp(x) 
a ' 



(4) 
(5) 



In the range < p < 2r the right spectrum of massless fermions is obtained [5]. 
For a quark of bare mass m the overlap-Dirac operator becomes 

(l - —am^)D N + m Q . (6) 



Since additive mass renormalization is forbidden, one avoids altogether a 
source of systematic errors that is always present with Wilson fermions [6]. 

Liischer has shown [7] that a fermion obeying the Ginsparg- Wilson relation 
possesses an exact chiral symmetry of the general form 

^ = e -75(1-- aL>W 5$ = e-^(l-^—^aD)-f 5 . (7) 
v p ' v p ' 
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A fully legitimate form of global chiral symmetry is thus preserved for finite 
lattice spacing. The expected value of the global anomaly [7] (which comes 
from the non-invariance of the fermionic integration measure under the trans- 
formations above), the Atiyah-Singer index theorem on the lattice [2] and the 
chiral Ward identities are then fully attained before taking the continuum 
limit. Thus, the central point in this framework is not to insist on the lattice 
form of the canonical chiral transformations (in fact actions satisfying Eq. (1) 
are not chirally invariant in the canonical sense), but rather be satisfied with 
a modified form of the chiral symmetry. The Nielsen- Ninomyia theorem [8] 
is then circumvented altogether, and it is possible to define chiral fermions 
without doublers or breaking of flavor symmetry or other unpleasant draw- 
backs. The chiral transformations (7) depend also on the interaction, but this 
does not forbid the non-perturbative construction of chiral gauge theories on 
the lattice with exact gauge invariance [9,10]. For recent reviews on the in- 
teresting properties of fermions satisfying the Ginsparg- Wilson relation, see 
Refs. [5,10,11]. 

We will be interested in the following in performing chiral-invariant compu- 
tations by using the overlap-Dirac operator (2). This has been proven to be 
local [12], and although simulations with the Neuberger operator look com- 
putationally very demanding, progress is under way [12,13]. There has been 
activity also on the analytic side, and some 1-loop calculations with overlap 
fermions have been already carried out [14-17]. The most recent and advanced 
calculations have featured the relation between the A parameter in the lattice 
scheme defined by the overlap operator and in the MS scheme [16], and the 
renormalization factors of the quark bilinears tpTip [17]. 

In this paper we present the calculation in lattice QCD of the renormalization 
factors of a few operators which measure the lowest two moments of non- 
singlet quark distributions. The generic operators which measure their n-th 
moments are ipj^D^ . . . D^ip and rfj^D^ . . . D^ n ip for unpolarized and 
polarized Structure Functions respectively. Contrary to what happens with 
Wilson fermions, in the present case, thanks to the exact chiral symmetry that 
we maintain, the renormalization constants for every pair of corresponding 
unpolarized and polarized operators which differ by a 75 matrix are the same. 

This paper is organized as follows: in Sect. 2 we introduce the various oper- 
ators that we have studied, in Sect. 3 their renormalization on the lattice is 
discussed, and in Sect. 4 some details about the perturbative calculations are 
given. Finally, in Sect. 5 we present our results. In the Appendices one can 
find the Feynman rules that we have used as well as some analytic results. 



The locality is not meant here to be strict locality, but in the larger sense that 
the strength of the interaction decays exponentially with the distance. 
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2 Moments of Structure Functions 

The hadronic physics in Deep Inelastic Scattering is contained in the matrix 
elements of the T-product of the hadronic currents 

-ijd 4 xe^{p\T(J,(x)J u (0))\p), (8) 

q being the momentum transfer and p the target momentum. An operator 
product expansion on the light cone of the kind 

J(x) J(0) ~ £ C^\x 2 ) x» ■ ■■x^% n (0) (9) 

n,i,l 

describes the physics in the Bjorken limit, in which scaling is reached and the 
Structure Functions depend only on the Bjorken variable xb = — q 2 /(2p ■ q). 
The moments of the Structure Functions, which measure the moments of quark 
distributions, are directly connected to the forward matrix elements of the 
local operators appearing in the light-cone expansion. Since the divergence of 
the Wilson coefficients is governed by the twist (dimension minus spin) of the 
corresponding operators, infinite towers of operators of increasing dimension 
and spin but with the same twist (dimension minus spin) appear at a given 
order in l/q 2 . The dominant contribution is given by twist two, which in the 
flavor non-singlet case means the symmetric traceless operators 

A a 

^..^^^A'-^my^ ( 10 ) 

O^^^i^D^-.-D^ij, (11) 

where A a are flavor matrices (which in the following will be omitted and im- 
plicitly understood). These operators measure the moment of unpolarized and 
polarized Structure Functions respectively, that is the moments (x B ) of quark 
momentum distributions and the moments ((Axb)™) of quark helicity distri- 
butions, at leading twist. Higher- twist operators give the sub-dominant con- 
tributions, and in particular twist-4 operators (which include also 4-fermion 
operators) measure the l/q 2 corrections to these moments [18]. 

Since we are considering flavor non-singlet operators, there can be no mixing 
with operators like Tr £ p F^ p D^ ■ ■ ■ F piln and Tr^F^D^ • • ■ F piln which 
measure the gluon distributions. Mixing with these operators is prohibited 
even when one considers unquenched calculations. Singlet quark operators 
instead do mix with gluon operators. For more detailed discussions of Deep 
Inelastic Scattering on the lattice, see Refs. [19-23]. 
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The operators in the light-cone expansion (9) of which we compute the renor- 
malization factors are the following: 



V2 , d = 'ih {1 D A yi/> (12) 
1 3 - 

0„ 2 , e = ^74-0^ - o S^TiW, (13) 

6 1=1 

which measure the first moment of quark momentum distributions, 



O a2 ,d = '4>l{il5D 4} iP (14) 
1 3 - 

O a2 ,e = 1pl4l5D 4 lfj - - ^ ^TiTsA^, (15) 

6 i=l 

which measure the first moment of quark helicity distributions, 



V3 ,d = i>l{ADiD 2 }ii (16) 

1 3 - 

V3 , e = "^E ^7{4AA}^, (17) 

which measure the second moment of quark momentum distributions, and 



O a3 ,d = ^7{47sA A}V> (18) 

1 3 - 

0« 3 ,e = ^7{47sA A}^ - - ^7{475 A A}^, (19) 

A i=2 

which measure the second moment of quark helicity distributions. Symmetriza- 
tion in all Lorentz indices is understood. We use D=D — D, with the following 
lattice discretizations: 



A i/)(x) — — U(x, fi)ip(x + afi) — U\x — afi, fi)ip(x — afi) (20) 

Zi(X 

ip(x) A = — + Q>p)U^(x, fi) — ip(x — afi)U(x — afi, //) . (21) 



We have chosen the Lorentz indices of each operator appearing in the con- 
tinuum expansion in two different ways, so that they fall in two different 
representations of the hypercubic group (the symmetry group of the lattice, 
remnant of the Lorentz symmetry), and on the lattice they will then renormal- 
ize in a different way. The representations where the indices are all different 
from each other are least likely to mix with other operators, however one 
needs more components of the hadron momentum different from zero when 
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simulating the corresponding matrix elements, and this can lead to more lat- 
tice artifacts. Sometimes one of the representations can be easier to handle in 
practice, depending on the trade-off between the amount of mixing and the 
number of non-zero momentum components. For example, in the case of the 
second moment the operator (16) is multiplicatively renormalized, but it can 
be more advantageous to use the other representation (operator (17)) when it 
is important that fewer components of the momenta are different from zero, 
although in this case one has then to deal with a mixing. 

The operator (17) (together with the corresponding polarized (19)) is the only 
one in the list above that is not multiplicatively renormalized on the lattice. 
What happens is that that the two operators 



1 3 

A = i)^D x D x i) V^A AV>, (22) 

B = ^ 7 i AAV> + i/S 7 i A AV> - ^ ]T ^D 4 D^ - ^ $H D l D A ^, (23) 

1 i=2 1 i=2 

which are not symmetrized in their indices, do not receive the same 1-loop 
corrections on the lattice [20], and so the operator (17), 

V3 , e = \(0 A + B ), (24) 



does not go into itself under lattice renormalization, because the symmetric 
combination is lost. In the polarized case, an analogous situation occurs for 



A s = ^74 7 5 A A^ - - ^7475 A Diip, (25) 

1 i=2 

B s = ^^DiD^ + ^7175-^1-04^ (26) 

1 3 - 1 3 - 

2 i=2 2 i=2 

A detailed discussion of these effects can be found in Sect. 5. 

We mention that for the second moment there exists also a third independent 
lattice representation, in which all indices are equal. However, it mixes with a 
lower dimensional operator with a power-divergent coefficient (~ 1 /a 2 ), and 
so we do not consider it here. 
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3 Renormalization 



To relate the numbers obtained from Monte Carlo simulations to physical 
continuum quantities, a lattice renormalization of the relevant matrix elements 
is necessary. The connection is given by 

(or> = e («« - -£- 2 - ) ) • > (27) 

where 

(0 co„t,la t) = £ ^. + jL^contaat^J . ^tree) (2g) 

are the continuum and lattice 1-loop expressions respectively, and the tree- 
level matrix element is the same in both cases. The differences Ai?^ = R 1 ^ 1 — 
Rf° nt enter then in the renormalization factors 

a 2 

Zij(afi, g ) = 5ij - Y^ARij(afx) (29) 

which connect the lattice to the continuum. These renormalization factors are 
independent of the state, depend only on the scale a\i and are gauge-invariant. 
Using them, a matrix element obtained from Monte Carlo simulations can be 
renormalized to a continuum scheme. 

Although there are in principle also non-perturbative methods with which one 
can determine the renormalization factors, perturbation theory is still impor- 
tant. In fact, it can happen that for some operators a window for the non- 
perturbative signal is difficult to obtain, and then perturbation theory remains 
the only possibility of computing the relevant renormalization factors. For the 
Neuberger operator, extracting the non-perturbative renormalization factors 
in addition to simulating the bare matrix elements could turn out to be com- 
putationally very expensive. In general, perturbative lattice renormalization is 
important as a hint and a guide for non-perturbative renormalization studies, 
and even more when mixings are present, which are generally more intricate 
than in the continuum case, and more transparent when looked at in pertur- 
bation theory, especially if some amounts of mixings are small. Perturbative 
renormalization can in any case be very useful in checking and understanding 
results obtained with non-perturbative methods. 

In the continuum we renormalize the operators in the MS scheme. As per- 
turbative renormalization condition on the lattice we impose that the 1-loop 
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amputated matrix elements at a certain reference scale [i are equal to the 
corresponding bare tree-level quantities. For lattice matrix elements of mul- 
tiplicatively renormalized operators computed between one-quark states this 
condition means 



(p\o l ^)\ P ) 



Zo(an,g ) ■ ZJ(an,g ) ■ (p\O m (a)\p) 



1— loop 



(p|O< >(a)|p) 



tree 



(30) 



where is the wave-function renormalization, computed from the quark self- 
energy. 

The 1-loop lattice matrix elements of the operators we consider here, which 
are at most logarithmically divergent, will have the form 



(p\0^(a)\p) 1 - loop = (p|0(°)(a)b> 



tree 



X 



(31) 



(l + ^2^(70 log aV + Vo + T + 2sfj, 



where Vq is the finite contribution of the vertex and sails diagrams (a, b and 
c in Fig. 1), T refers to the tadpole arising from the operator (d in Fig. 1), S 
is the finite contribution (proportional to ij>) of the quark self-energy of one 
leg, including the leg tadpole (e and g, or f and h, in Fig. 1), and Cp = ^^T 1 . 
The Z Q factor for an operator O will then be given by 



Z o (afi,g ) 



16tt 2 



C F I 70 log a V + B ), 



(32) 



with 

B = Vo + T + S. (33) 



We will call "proper" contributions the ones that exclude the self-energy dia- 
grams. They correspond to the diagrams a-d in Fig. 1. 

Although the theory defined by the overlap-Dirac operator (2) has an ex- 
act chiral symmetry and no lattice artifacts of order a are present in the 
action, matrix elements of operators still possess corrections of order a and 
therefore they need to be improved. The improved operator corresponding to 
V>7{^---^W is [24] 

^(l - ^aD N ) ■ ■ ■ D^y (l - ^aD N )^, (34) 
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(a) vertex (b) sail 




(c) sail (d) operator tadpole 




(e) leg self-energy (f) leg self-energy 

(excluding tadpoles) (excluding tadpoles) 




(g) leg tadpole (h) leg tadpole 



Fig. 1. The graphs contributing to the 1-loop renormalization factors of the matrix 
elements (p\0\p). The operator insertion is indicated by a circle. 
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and the big bonus with overlap fermions is that the renormalization factors 
for the improved operator and for the unimproved operator without the rota- 
tions (1 — -aDjv) are the same. What happens is that in 1-loop amplitudes a 
factor Djsf can combine with a quark propagator, but since it has an a in front 
and (contrary to the Wilson case) there is no l/a factor in the propagator, as 
additive mass renormalization is forbidden by chiral symmetry, the contribu- 
tion of Dn to the renormalization factors is zero [17]. Thus, we can simulate 
the improved operator (34) and renormalize it as if it were the unimproved 
operator, which saves a lot in the perturbative calculations. 

Chiral symmetry, in addition to avoiding the extrapolations to zero quark mass 
in the simulations, gives also useful constraints on the mixing patterns and on 
the values of the renormalization factors. In general, mixings that in the Wilson 
case were allowed by the breaking of chirality are now forbidden, although this 
has less relevance here than in other problems like in weak interactions [25]. 
An important consequence here is that, similarly to what happens in the case 
of the bilinears where we have Zg = Zp and Zy = Za [17], chiral symmetry 
forces pairs of corresponding unpolarized and polarized operators, like V2j a 
and O a2; d, to have the same renormalization constants 0. 



4 The computations 

The interaction vertices and the propagator of the overlap-Dirac operator are 
much more complicated than the ones in the Wilson formulation, and this 
causes the computations to be rather cumbersome (see Appendices). Com- 
puter programs need then to be introduced, and we have performed the cal- 
culations using an ensemble of routines written in the symbolic manipulation 
language FORM. These routines are an extension of the ones used to do cal- 
culations with the Wilson action in various occasions [23,26]. 

The outputs of the FORM codes, which correspond to the results of the an- 
alytic calculations, are fed to Fortran programs which perform the numerical 
integrations. To treat the divergent integrals we use the Kawai method [27], 
in which they are split in integrals independent of external momenta plus 
continuum integrals. Divergent lattice integrals which contain both gluon and 
overlap-quark propagators are reduced to divergent integrals containing only 
gluon propagators, which are computed in a exact way as in [28] , plus a finite 
rest. 



2 This is strictly true for improved operators like in Eq. (34), which transform like 
chiral multiplets, but since they have the same renormalization constants as the 
unimproved operators, these symmetry relations are preserved. 



10 



Many checks have been performed on our results. We did the computations 
in a general covariant gauge 0, and the cancellations of the gauge-dependent 
parts between the continuum scheme and the lattice results is a strong check 
on the good behavior of the FORM codes, as well as of the integration routines. 

We have used the codes to compute also the renormalization constants for 
all quark bilinears, and we agree with the results given in Feynman-gauge by 
Alexandrou et al. [17], as well as with their self-energy. In the case of the leg 
tadpole and of the vertex graph of the scalar current ipip we have also made 
the computations entirely by hand (although only in the Feynman gauge), 
and successfully checked them against the FORM output expressions. These 
analytic results are reported in Appendix B. We tried to do a calculation by 
hand for a first moment operator and for the rest of the leg self-energy, but they 
involved a huge amount of manipulations. We rely also on the vast amount 
of results which have been produced in different occasions using these FORM 
codes limited to the Wilson formulation [23,26], and these results have by now 
a certain number of cross-checks. With overlap fermions we need to introduce 
a new propagator and new vertices, but the operators are unaltered, and the 
same Wilson expressions for their expansion in a can be used here. Many of 
the routines (for example the gamma algebra reduction) are also exactly the 
same as in the Wilson case. 

We checked in all cases that the polarized operators have the same Zs as 
the corresponding unpolarized operators. This is a rather strong check, as 
the analytic results are widely different for the two cases, and it is only after 
performing the numerical integrations and seeing that the difference of the 
two results is much smaller than the precision of our integrals that we able to 
say that the renormalization factors are indeed equal. Finally, since our codes 
are able to do calculations both in Dimensional Regularization and using a 
mass regularization, we also checked that using two different intermediate 
regularizations while treating the divergences with the Kawai method [27] we 
get the same results. 

We obtain the Bq constants with 5 digits after the decimal points, by doing the 
numerical integrations on a 40 4 regular grid and using the method proposed 
by Liischer and Weisz in [29] to accelerate convergence. We have checked that 
these digits do not vary when a 60 4 integration grid is used. In any case, to get 
more significant digits, grids of 80 4 or 100 4 points would still not be enough, 
and they would require an enormous computational effort, as a typical FORM 
output for the second moment operators already contains thousands of terms. 

3 The gluon propagator that we use is 

G - ( ^lE^¥( v - (1 - a) lM)' (35) 
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Comparing our results for the self-energy and the bilinears with the numbers 
in Ref. [17], we can see that a few times a difference of one unit on the fifth 
digit can be noticed, and this agrees with our estimate of errors. In that paper 
the results are given in terms of the constants b = and this explains the 
presence of two more significant digits after the decimal point. 



5 Results 



We give in this Section the renormalization factors of the operators consid- 
ered in this work. That these factors are identical for pairs of corresponding 
unpolarized and polarized operators has been explicitly verified for all oper- 
ators. We give then the numerical results only for the unpolarized operators, 
for r = 1. 

We first consider the 1-loop contributions of the proper diagrams (a-d in 
Fig. 1), which for the operators that do not mix are 

(J + (l-a))logaV 

+Vv 2 a t - (1 - «) 6.850272 + T V2 , d 
(| + (l-a))logaV (36) 

+V V2 a = 1 - (1 - a) 6.850272 + T V2 , e 
(y + (!-«)) log a 2 / 

+V V3 a f -(1- a) 7.369693 + T V3>d 

The contributions of the sails and vertices, Vo, have been for convenience 
separated in the Feynman gauge results (V^f 1 , V V2 7 X an d K^T 1 )) tabulated 
in Table 1 for various values of the parameter p, and the remaining parts 
proportional to (1 — a), which are instead independent of p. The analytic 
expressions of the latter are very complicated functions of p containing thou- 
sands of terms, and the numerical cancellation of this dependence is a rather 
non-trivial check of our computations. These numbers are also independent of 
the lattice representation, and furthermore they have the same value as with 
the Wilson action [26], so they seem to be to a certain extent independent of 
the particular fermion action chosen. 



^proper 9o ^ t 



V2,< 



167T 2 ~ F 



^proper 9p 



V2,e 



16tt 2 



^otree 



^proper 9p ^ ^ 



V3,' 



167T 2 



U v 3 ,d- 



12 



p 


V a 7 l = V "r 1 


V a=1 = V a=1 


V a j l = V a j l 


0.2 


-3.77889 


-3.50634 


-7.57984 


0.3 


-3.71038 


-3.41915 


-7.45926 


0.4 


-3.65209 


-3.34178 


-7.35607 


0.5 


-3.60037 


-3.27058 


-7.26432 


0.6 


-3.55327 


-3.20357 


-7.18078 


0.7 


-3.50959 


-3.13952 


-7.10341 


0.8 


-3.46851 


-3.07762 


-7.03088 


0.9 


-3.42946 


-3.01727 


-6.96221 


1.0 


-3.39203 


-2.95803 


-6.89668 


1.1 


-3.35589 


-2.89955 


-6.83375 


1.2 


-3.32077 


-2.84155 


-6.77296 


1.3 


-3.28647 


-2.78380 


-6.71398 


1.4 


-3.25281 


-2.72610 


-6.65650 


1.5 


-3.21964 


-2.66827 


-6.60029 


1.6 


-3.18685 


-2.61015 


-6.54513 


1.7 


-3.15430 


-2.55162 


-6.49084 


1.8 


-3.12190 


-2.49255 


-6.43729 



Table 1 

The Feynman-gauge constants V a=1 for the multiplicatively renormalized opera- 
tors. 

The contributions Tq of the operator tadpoles (d in Fig. 1) are shown in Table 
2, which contains also their results for the other operator, V3fi . However, since 
this operator is not multiplicatively renormalized on the lattice, we consider 
in its place the two operators Oa and Ob introduced in Sect. 2. They mix 
with each other and we can write their renormalization as Q 



6 A = ZaaOa + ZabOb (37) 
Ob = ZbaOa + ZbbOb- 

In terms of the bare operators Oa and Ob, the renormalized operator V3tC 

4 In Refs. [21] another choice for the two operators that mix was made, and V3jC 
and an operator of mixed symmetry were considered instead of Oa and Ob- 
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operator 


operator tadpole 


V2 ,di O a , 2t d 
V3t di O aZ) d 


T V2 , d = T a24 = -167r 2 f(l - |(1 - a)) 
T V2>e = T a2 , e = -16vr 2 f^(l - |(1 - a)) 
T V3>d = T a3>d = 16tt 2 ( - Z + (1 - a) f ) 


VSt e, O a3 ^ e 


Taa = T A 5 A5 = 16^ 2 ( - 2Z + § + (1 - a)(z - §)) 
Tab = B 5 = 
Tba = T B 5a-> = 

T BB = ?W = 16vr 2 ( - Z + (1 - a) f ) 



Table 2 

The operator tadpoles for the various operators, where Zq = 0.154933390231. 
appearing in the DIS light-cone expansion will be 

d V3 , e = \(Oa + Ob) = \{Z A a + Z BA ) O a + \{Z AB + Z BB ) O b , (38) 



and since on the lattice Zaa + Zba 7^ Za B + Z BB , the result is that 0„ 3je is not 
multiplicatively renormalized. In fact, the 1-loop contributions coming from 
the proper diagrams are in this case 



^proper 

° A "16^ 



C r 



2„2 



a) ) log a p 



^proper 



+ v aa = - (1 - «) 8.685568 + T AA 
1 



Otree 
A 



log aY + V AB = L + -(l-a) + Tab 



^B 



16tt 2 



2 log aY + VgT 1 + -(l-a)+T BA 
(^ + (l-a))\ogay 



Otree 
A 



(39) 



+V BB =1 ~ (1 - a) 7.703026 + T BB 



A-)tree 
^B ) 



where the finite contributions of vertex and sails in the Feynman gauge are 
shown in Table 3. 

We stress again that this mixing is a pure lattice effect, deriving from the 
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p 


t/ a=l _ t/ a=l 


t/ a=l _ t/ a=l 
V S,4 ~~ V B 5 A 5 


t/ a=l _ t/ a=l 


t/ a=l _ y a=l 


0.2 


22.83163 


-36.93402 


-13.77026 


10.81988 


0.3 


13.03570 


-24.07048 


-8.80836 


4.50553 


0.4 


8.39467 


-17.86369 


-6.43263 


1.50207 


0.5 


5.76489 


-14.27357 


-5.06826 


-0.20460 


0.6 


4.11604 


-11.96971 


-4.19847 


-1.27644 


0.7 


3.01383 


-10.38860 


-3.60511 


-1.99319 


0.8 


2.24472 


-9.25173 


-3.18071 


-2.49274 


0.9 


1.69205 


-8.40616 


-2.86648 


-2.85062 


1.0 


1.28701 


-7.76116 


-2.62765 


-3.11148 


1.1 


0.98657 


-7.25963 


-2.44241 


-3.30329 


1.2 


0.76254 


-6.86397 


-2.29645 


-3.44445 


1.3 


0.59577 


-6.54843 


-2.18002 


-3.54747 


1.4 


0.47280 


-6.29483 


-2.08624 


-3.62114 


1.5 


0.38399 


-6.08998 


-2.01015 


-3.67183 


1.6 


0.32223 


-5.92413 


-1.94810 


-3.70420 


1.7 


0.28218 


-5.78989 


-1.89733 


-3.72177 


1.8 


0.25980 


-5.68156 


-1.85572 


-3.72723 



Table 3 

The Feynman-gauge constants V a=1 for the operators Oa and Ob related to O. 



breaking of the Lorentz group to the hypercubic group. In the continuum we 
do have Zaa + Zba = Zab + Zbb, and thus the operator 0„ 3ie is multiplica- 
tively renormalized, and has the same Z as the operator with different indices 
V3t d, as they both belong to the same representation of the Lorentz group. 
We remind also that in the polarized case one encounters exactly the same 
situation with the same numerical results, i.e. Zaa — Z a ^a 5 , Zab — Za$b 5 
etc., as we have explicitly verified. 

To complete the computation of the renormalization factors, we have now 
to add to the proper diagrams the 1-loop amplitudes of the self-energy and 
tadpole of one leg which are proportional to ift, 



a log a V + S a=L + (1 - a) 4.792010 



(40) 
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p 


S 


P 


S 


P 


S 


0.2 


-240.59963 


0.8 


-49.99380 


1.4 


-23.76598 


0.3 


-155.41069 


0.9 


-43.10649 


1.5 


-21.50111 


0.4 


-112.98999 


1.0 


-37.63063 


1.6 


-19.53531 


0.5 


-87.65282 


1.1 


-33.17935 


1.7 


-17.81537 


0.6 


-70.84428 


1.2 


-29.49505 


1.8 


-16.29999 


0.7 


-58.90122 


1.3 


-26.39961 


1.9 


-14.95658 



Table 4 

1-loop results for the quark self-energy (including the leg tadpole). 



where the Feynman-gauge finite results S a=1 are given in Table 4. Putting 
everything together, we get the expressions of the renormalized operators on 
the lattice for overlap fermions, which for p = 1 are: 



Ov 2 ,d ■ 



1 - 
1 - 



j^C F ^ log a V - 53.25571 + (1 - a) 
-^- 2 C F f| log a 2 /i 2 - 52.82171 + (1 - a) 



O 



tree 
v 2 ,d 



s~\tree 

v 2 ,e 



1 - ^C F (j log a 2 /i 2 - 68.99341 + § (1 - a)) 

l07T Z V D 2 



.otrcc 



0\ rcc (41) 



+- 



1 - rf^CV f ^ log a 2 /i 2 - 67.83586 + ^ (1 - a) 
\Q-n z \ 6 2 



Otree 



We give here for comparison the numerical values obtained with the usual 
Wilson action, without any improvement. We have computed them again and 
checked them with the results in the literature [19-23,26]. In this case however 
chiral invariance is broken and the renormalization factors for the unpolarized 
operators 



,0 Wilson 
U V2,d — 


1 - 


9l 
16tt 2 


/^Wilson 


1 - 


9l 
16tt 2 



-C f - log a 2 // - 3.16486 + (1 - a) 
V 3 



,2, .2 



U v 2 ,d 



O 



tree 
V2,e 
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/^Wilson 
U V 3 ,d 



o 



Wilson . 

v 3 ,e 



1 - 
1 



+ - 



1 



JjLc, (f log aV - 19.00763 + § (1 - a)) 

- i^Cr^ log aW - 2178271 + | (1 - a) 

9o q (25 n 2 2 _ „ .„„.„ 3 



C^ ree (42) 



1 



. ( .__ 2 ^ F . — log a 2 p 2 - 18.46640 + -(!-«) 



differ from the ones of the polarized operators 



A-itree 



/^Wilson 

U a 2 ,d — 

/^Wilson 

/^Wilson 

U a 3 ,d - 

/^Wilson ^_ 



1 _ Jo 



j^C F log a 2 fi 2 - 4.09933 + (1 - a) 
^C F ^ log a 2 fi 2 - 4.27705 + (1 - a) 



.otree 



^jtree 

Q2,e 



g2 ° -C F ( ^ log aV " 19.56159 + 1(1 -a) 
6 2 

1 -ifMir logaV - 



16tt 2 



.otree 
U a 3 ,d 



1 

+ 3 



25 3 
— log a 2 ^ 2 - 22.39940 + - (1 - a) 
6 2 

19.25837 +|(l-a) 



1 -ilHf logaV - 



O^r (43) 



Otree 
£5 • 



What one can notice is that the values of the renormalization factors for 
overlap fermions are in general much larger than for Wilson fermions. This 
seems to be true for most values of p, and is to trace largely to the quark 
self-energy and to the operators tadpoles. It can be observed that when in 
the Wilson formulation one improves the theory by adding the clover term to 
the action (with c sw = 1) and by canceling 0(a) effects on the operators at 
tree level, the renormalization factors become also somewhat larger, as in the 
examples below [19,20]: 



^Wilson, imp 
U V2,d 



o 



Wilson, imp 

V3,d 



1 - 



-|Lc F (|logaV- 15.816 + (1- a) 
9o q f 25 , ^ j ..2 - 3 



^trcc 
U V 2 ,d 



16tt 2 



/^Wilson, imp 

■ — g 



+ - 



1 - 



9l 



167T 2 
1 - 



6 log ay - 29.815 + -(1- a) J 

/25 3 
C F — log a 2 /i 2 - 39.192 + - (1 - a) 
\ 6 2 



.otree 
U v 3 ,d 



9o q f 25 , ^ n2 , .2 — - 3 



16tt 2 



log a 2 u 2 - 22.141 + -(1- a) 

6 2 v ' 



0^44) 



/Qtrcc 



One could then speculate whether the large renormalization factors that we 
have obtained for overlap fermions are related to the fact that the Neuberger 
action is improved and the Zs correspond to the ones of improved operators. 
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A calculation with Wilson fermions with which to compare in this sense our 
overlap results would be one in which the operators are improved beyond tree 
level, but although we know in some cases the contributions of the operator 
counterterms that are needed for the full improvement, unfortunately not all 
their coefficients are yet determined [26]. 

We give finally the 1-loop results for the various matrix elements in the MS 
scheme, so that one can complete the connection with the continuum as in 
Eq. (27). In the continuum there is only one renormalization constant for all 
possible operators (unpolarized and polarized) measuring the first moment, 
and only one for all possible operators measuring the second moment: 



/OMS _ 

V2 


1 - 


9l 
16tt 2 


nMS_ 
w v 3 


1 - 


9 2 o 
16tt 2 



_ /8. v 40 M 



^otree 

v 2 



(45) 



25 
6 



log 



P 
J 2 



67 3 M 

y + 2 (1 



a) 



jQtree 
W V3 • 



The connection of overlap lattice fermions with the continuum MS is then 
given by the gauge-invariant factors 



/=jMS 

U V 2 ,d 



o 



MS . 

V2,e ' 



/=jMS 
U v 3 ,d 



O 



MS . 

v 3 ,e ' 



C F [ ^ log aV -48.81127 



1 - 
1 - 



16tt 2 1 \3 
I |^ci^logaV-48.37727 N 



V2,e 



^cJ | log aV" 61.54897) 



^lat 



1 

3 



1 - 



go 
16tt 2 



C F f — log a 2 fi 2 - 65.85333 ) 



6 



l 

+ 3 



^C^l log «V 2 - 60.39142 j 



,olat 



(46) 



It looks increasingly difficult to go to higher moments, as the number of terms 
that are present in the FORM outputs and that need to be numerically inte- 
grated becomes very large. This has at the moment limited our calculations 
to second moment operators, but we hope to be able to compute the renor- 
malization of third-moment operators in the near future. 

Finally, we remark that all the renormalization constants presented here can 
be considered as computed in the unquenched case, because we limit ourselves 
to 1-loop computations with non-singlet quark operators, where internal quark 
loops never have the chance to come to play. 
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A Feynman rules for the Neuberger operator 



We give here the Feynman rules which are needed to perform 1-loop calcula- 
tions using the overlap-Dirac operator. An explicit derivations of these rules 
is given in [14,15]. 



The quark propagator is 



b{k) ~ 2p(u;(k) + b(k)) + 2~p' {AA) 



where 



2 



u(k) — —a sin 2 ak u + f 2r sin 2 — — — p 
W = i(2rEsm 2 ^-p). (A.2) 



The vertices needed for 1-loop calculations can be entirely given in terms of 
the vertices of the Wilson action 



W lll {p u p 2 ) = ~g \t^ ll cos + rsm I (A.3) 

W2^Pi,P2) = -^agA -i^sm ^ + rcos (A.4) 

(where p\ and p 2 are the quark momenta flowing in and out of the vertices) 
and of the quantity 

Xo{p) = - U^7M sm «PM + 2r H sm2 -TT ~ p]- ( A - 5 ) 
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The quark-quark-gluon vertex in the overlap theory has the expression 



Vi li (pi,p 2 )=p 



1 

u(pi) +uj(p 2 ) 



x 



and the quark-quark-gluon-gluon vertex is 



Xo( P2 )Wl(p 1 , P2 )X (p 1 ) 



(A.6) 
, (A.7) 



V 2 ^(pi,p 2 ) =6^p- 



-x 



W 2 ^ Pl ,p 2 ) 1 X Q (p 2 )Wl{p l ,p 2 )X (p 1 ) 

UJ{Pl)UJ{p 2 ) 

11 1 1 

P 

2 uo{pi) + u{p2) v(pi) + v(k) u;(k) + uj(p 2 ) 

X ( P 2)Wl(p 2 ,k)W lu (k,p 1 ) + W l ,(p 2 ,k)X} ) (k)W lu (k,p 1 ) 

+W^(p 2 ,k)wL(k,p l )X (p 1 ) 
u(pi) + u(k) + Ld(p 2 ) 



(A.8) 



u(pi)u(k)u(p 2 ) 



X (p 2 )Wl(p 2 ,k)X (k)Wl(k, Pl )X ( Pl ) 



B Some analytic results 



We give here the analytic results for the leg tadpole and for the vertex of the 
scalar current, in the Feynman gauge for r = 1. In order to be able to write 
them in a compact form it is convenient to introduce the abbreviations 



k\ k\ 
M A = cos— , A" A = sin— , s\ = sin k\, s 2 = J2 s li 



5 = 5(A0 = 2£sin 2 ^-p, 



D = 2p(uj(k) 



A 



u 2 {k) 



2 k\ 



1 sin 

Pa 2 



' V A \ A 



sin 



2 



(B.l) 

(B.2) 
(B.3) 
.(B.4) 



The result for the 1-loop leg tadpole is then given by 
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X 



M 2 + iV 2 + (\ + ^^- s J + 6 (M 2 -iV 2 ) 
2 + ' 



E 

A 



+ 



+00 



G 2 (A;) 



- 6(M 2 - Ni) + S 2 



1 



2tt 4 v V(l + V^) 2 



x 



E 



- K N l + ( 1 + 7i ) K (6 + 2M l - M ^ 



(B.5) 



The first term comes from the part of the V2 vertex (A. 8) containing W% and 
W 2 \ and its value for p = 1 is 



Z "u I -1 



4 X 



9l 



p=i 



16tt 2 



36.69915, 



(B.6) 



while the 1-loop result for the whole leg tadpole is smaller: 



16tt 2 



23.35975. 



(B.7) 



Adding now the value — j^jp 14.27088 for the diagram e in Fig. 1, which is much 
harder to compute by hand and would be a very lengthy analytic expression 

2 

anyway, gives the result — yf^- 37.63063 for the complete self-energy in the 
Feynman gauge for p = 1. 



The results for the 1-loop vertex diagram of the scalar operator (the sails are 
not present in this case) can be written in the form 



9o 



s 

£2 



Ap\ 



X + — -Y 
pD 



(B.8) 



for small p, where 



*=E 



- (M 2 X - N 2 
1 



26(M 2 + TV 2 ) 



p 2 A 



pVA 

s 2 -b 2 ){Ml -Nl) + 2bs 2 



21 



2s 2 (M 2 x + N 2 X ) 



p 2 A 



s 2 [ -2b(M 2 -N 2 ) + s 2 -b 2 



(B.9) 



We have checked that all these results obtained by hand correspond with the 
outputs of the FORM programs. 
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